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Based on a system-reservoir model, where the system is nonlinearly coupled to a heat bath and the 
heat bath is modulated by an external stationary Gaussian noise, we derive the generalized Langevin 
equation with space dependent friction and multiplicative noise and construct the corresponding 
Fokker-Planck equation, valid for short correlation time, with space dependent diffusion coefficient 
to study the escape rate from a metastable state in the moderate to large damping regime. By 
considering the dynamics in a model cubic potential we analyze the result numerically which are in 
good agreement with the theoretical prediction. It has been shown numerically that the enhancement 
of rate is possible by properly tuning the correlation time of the external noise. 
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I. INTRODUCTION 

Barrier crossing phenomena are ubiquitous and are 
a central issue in many areas of natural science 0, 
Since the reaction coordinate describing the transition 
across the barrier typically interact with a large num- 
ber of microscopic degrees of freedom, a useful the- 
ory to start with is the construction of a Hamiltonian 
which incorporates the effects of a heat bath environ- 
ment. The classical treatment of thermally activated 
barrier crossing description is naturally accounted for by 
the generalized Langevin equation (may be nonlinear) 
and by the corresponding Fokker-Planck equation. Fol- 
lowing Kramers 0, 01 transition rate can then be cal- 
culated from a nonequilibrium steady state solution of 
the Fokker-Planck equation describing a constant flux 
across the potential barrier. Over several decades the 
Kramers' theory and many of its variants has served as 
standard paradigms in various problems of physical and 
chemical kinetics to understand the rate in multidimen- 
sional systems in the overdamped and underdampcd lim- 
its effects of anharmonicities Q, rate enhancement 
by parametric fluctuations [1] , the role of non-Gaussian 
white noise 0, , role of a relaxing bath [lo| , quantum 
and semiclassical corrections to classical rate processes 
[O, and related similar aspects. 

The common feature of overwhelming majority of the 
aforesaid treatments is that the system is thermodynam- 
ically closed, which means that the noise of the medium 
is of internal origin so that the dissipation and fluctu- 
ations gets balanced through the fluctuation-dissipation 
relation. However, in a number of situations the system is 
thermodynamically open, i.e., when the system is driven 
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by an external noise which is independent of system's 
characteristic damping jT3|. The main feature of the 
dynamics in this case is the absence of any fluctuation- 
dissipation relation. While in the former case a zero cur- 
rent steady state situation is characterized by an equi- 
librium Boltzmann distribution, the corresponding situ- 
ation in the later case is defined only by a steady state 
condition, if attainable. 

A common approach to study the nonlinear, nonequi- 
librium systems involves a description in terms of non- 
linear stochastic differential equation [Tsj . From a mi- 
croscopic point of view, the system-reservoir Hamilto- 
nian description suggests that the coupling of the sys- 
tem and the reservoir coordinates determines both the 
noise and the dissipative terms in the Langevin equation 
describing the dynamics of the system. If the system- 
bath interaction is linear in the bath coordinates but 
arbitrary in the system coordinate, the corresponding 
generalized Langevin equation incorporates a multiplica- 
tive noise and consequently a nonlinear dissipative term 
arises due to the nonlinear system-bath interaction p^ . 
A canonical distribution for initial conditions of the bath 
variables yields a zero average for the fluctuations and 
the fluctuation-dissipation relation is again maintained 
[l6| . However, when the reservoir is modulated by an 
external noise, it is likely that it induces fluctuations in 
the polarization of the reservoir • Since the fluctua- 
tions of the reservoir are crucially dependent on the re- 
sponse function, one can envisage a connection between 
the dissipation of the system and the response function 
of the reservoir due to the external noise from a micro- 
scopic stand point [l^- A direct driving of the system 
usually breaks the fluctuation-dissipation relation and 
can generate a biased directed motion that are seen in 
ratchet and molecular motors 18]. On the other hand the 
bath modulation by an external noise agency maintains a 
thermodynamic consistency relation, an analogue of the 
fluctuation-dissipation relation of the closed system, as a 
result of which the well known Kramers' turnover feature 
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can be restored [17|. 

While the nonequihbrium, nonthermal systems have 
also been investigated phenomenologically by a number 
of works in several contexts [l9|, these treatment are 
mainly concerned with direct driving of the system by an 
external noise or a time dependent field. In the present 
paper we consider a system-reservoir model where the 
reservoir is modulated by an external noise and the sys- 
tem is nonlinearly coupled to a heat bath thereby re- 
sulting a nonlinear, multiplicative generalized Langevin 
equation. Our object here is to explore the role of reser- 
voir response on the system's dynamics and to calculate 
the generalized escape rate from a metastable state for 
a nonequilibrium open system in the presence of multi- 
plicative noise. 

A number of different situations depicting the modula- 
tion of the heat bath may be physically relevant. As for 
example, one may experimentally study the reversible 
isomerization of czs-butene to irans-butene, cis-butene 
^ trans -hutene. In terms of the reaction rate theory 
both isomers represent the two stable local minima of 
the potential energy landscape and are separated by the 
activation energy which one stable configuration {cis or 
trans) needs to overcome to be converted into its isomeric 
form. To observe the effect of external stochastic modu- 
lation one can carry out the experiment in a photochem- 
ically active solvent (the heat bath) where the solvent is 
under the influence of external monochromatic light with 
fluctuating intensity of a wavelength which is absorbed 
solely by the solvent molecules. As a result of it the mod- 
ulated solvent heats up due to the conversion of light en- 
ergy into heat energy by radiationless relaxation process 
and a effective temperature like quantity develops due 
to the constant input of energy. Since the fluctuations 
in the light intensity results in the polarization of the 
solvent molecules, the effective reaction field around the 
reactants gets modified. Provided the required station- 
arity of this nonequilibrium open system is maintained 
the dynamics of the barrier crossing becomes amenable 
to the present theoretical analysis that follows. 

Though the dynamics of a Brownian particle in a uni- 
form solvent is well-known, it is not very clear when the 
response of the solvent becomes time dependent, as in a 
liquid crystal when projected onto an anisotropic stochas- 
tic equation of motion or in the diffusion and reaction in 
supercritical liquids and growth in living polymerization 
[20I |2]| . Also space-dependent friction is realized from 
the presence of a stochastic potential in the Langevin 
equation (2^ . An exact Fokker-Planck equation for time 
and space dependent friction was derived by PoUak et al 
[2^ 1. Along with these formal developments, the theories 
of multiplicative noise have found in wide applications 
in several areas, e.g., activated rate processes, stochas- 
tic resonance, laser and optics, noise induced transport 
In passing we mention that the escape rate for 
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a space dependent friction is just not a theoretical issue 
but has been a subject of experimental investigation over 
the last two decades 1251. 



The organization of the paper is as follows: In sec- 
tion II, starting from a Hamiltonian description of a sys- 
tem nonlinearly coupled to a harmonic reservoir which 
is modulated externally by a Gaussian noise, we have 
derived the generalized Langevin equation with an ef- 
fective Gaussian noise £^{t). The statistical properties of 
^(t) has also been explored. In section III we have con- 
structed the corresponding Fokker-Planck equation, valid 
for small correlation time, and have derived the general- 
ized Kramers' rate for moderate to large friction. In sec- 
tion IV, a specific example has been carried out. Both 
the numerical and analytical results have been analyzed 
in section V. The paper has been concluded in section 
VI. 



II. THE MODEL: HEAT BATH MODULATED 
BY EXTERNAL NOISE 

We consider a classical particle of mass M nonlinearly 
coupled to a heat bath consisting of N harmonic oscilla- 
tors driven by an external noise. The total Hamiltonian 
of such a composite system can be written as [26|, l23| 



H = 



P 



1 



N 



+Hint 



p 



m 



Cjg{x)y 



(2.1) 



where x and p are the coordinate and the momentum of 
the system particle, respectively, and V{x) is the poten- 
tial energy of the system, {qj , pj ) are the variables for the 
jth bath oscillator having frequency coj and mass rrij . Cj 
is the coupling constant for the system-bath interaction 
and g{x) is some analytic function of system coordinate. 
Hint is the interaction term between the heat bath and 
the external noise e{t), with the following form: 



Hi, 



N 



qje{t). 



(2.2) 



The type of interaction we have considered between the 
heat bath and the external noise. Hint is commonly 
known as the dipole interaction j28| . In equation (|2.2p . 



Kj denotes the strength of the interaction. We consider 
e(t) to be a stationary, Gaussian noise process with zero 
mean and arbitrary correlation function 



(e(i))e = 0, (e(i)e(i'))e = 2Z?*(t-i') 



(2.3) 



where D is the external noise strength, ^'(i — i') is the 
external noise kernel and (...)e implies averaging over the 
external noise processes. 

Eliminating the bath degrees of freedom in the usual 
way [16[ (and setting M and nij = 1) we obtain the 
generalized Langevin equation 
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dV{x) dg{x) 
dx dx 



(2.4) 



where 



AT 

7W ^^c^uj] cos ivjt 



(2.5) 



and /(t) is the thermal fluctuation generated due to the 
system-reservoir interaction and is given by 



N 



/(*) = X! [^J^H^jW ~ "^Jf (2^(0))} COS cjji 

«7(0) . 
H smud . 



(2.6) 



In equation (|2.4p . 7r(t) is the fluctuating force generated 
due to the external stochastic driving e{t) and is given 
by 



{{m^{t'))) - kBT'i{t-t') + 2D f dt" f dt'" 

Jo Jo 

Xip{t - t")ip{t' - i"')*(t" - t'") (2.9) 

along with {{£,{t))) = 0, where ((...)) means we have taken 
two averages independently. It should be noted that the 
above equation (|2.9p is not a fluctuation-dissipation rela- 
tion due to the appearance of the external noise intensity. 
Rather it serves as a thermodynamic consistency condi- 
tion. The statistical properties of n(t) are determined 
by the normal mode densities of the bath frequencies, 
the coupling of the system with the bath, the coupling 
of the bath with the external noise and on the statisti- 
cal properties of the external noise itself. Equation (|2.7p 
is reminiscent of the familiar linear relation between the 
polarization and external field, where 7r(i) and e{t) play 
the role of the former and later, respectively. ip{t) then 
can be interpreted as a response function of the reservoir 
due to the external noise e{t). The structure of 7r(t) sug- 
gests that this forcing function, although obtained from 
an external agency, is different from a direct driving force 
acting on the system. 

To obtain a finite result in the continuum limit, the 
coupling function Ci = c(w) and Ki — k{u!) are chosen 17 1 
as c{(jj) = co/wy^ and k{uj) = n^uj^fT'c. Consequently 
7(t) and </3(t) reduces to the following forms: 



where 



■K{t) = - / dtV(i - i )e(i ) 

^0 



(2.7) 
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dLDSl{^J) cos 



(2.10) 
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CjUjKj sin LOjt. 



(2.8) 



The form of equation (|2.4p indicates that the system is 
driven by two forcing terms f(t) and 7r(t), both are multi- 
plicative by a function of system variable dg{x) / dx. Thus 
we have obtained a generalized Langevin equation with 
multiplicative noise. To define the statistical properties 
of f(t), we assume that the initial distribution is one in 
which the bath is equilibrated at i = in the presence 
of the system but in the absence of the external noise 
agency e{t) such that 



{fit}) = and {fit) fit')) = keTjit-t'). 

Now at i = 0+, the external noise agency is switched on 
and the bath is modulated by e(t). The system dynam- 
ics is governed by equation l|2.4p . where apart from the 
internal noise fit) another fluctuating force nit) appears 
that depends on the external noise e(t). So we define 
an effective noise ^(i)[= fit) + 7r(t)] whose correlation is 
given by [l3| 



(fit) = cqKq / dLu!^iuj)LLj sin ujt 



(2.11) 



where co and kq are constants and I/tc is the cutoff fre- 
quency of the oscillator. Tc may be characterized as the 
correlation time of the bath [l^. For Tc ^ we ob- 
tain a delta-correlated noise process, i^(ijj) is the den- 
sity of modes of the heat bath which is assumed to be 
Lorcntzian: 



^iiv) = -- ^ 



TTTciuj'^ +Tc ^) 



(2.12) 



This assumption resembles broadly the behavior of the 
hydrodynamical modes in a macroscopic system [20| . 
With these forms of ^(w), c(a;) and niuj) we have the 
expression for ipit) and "/{t) as 



'Pit) 
lit) 



Tc 
^2 



exp(-i/Tc) 



^exp(-Vrc). 



(2.13a) 
(2.13b) 
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From equations (|2.10p and (|2.1ip we obtain 



dt 



Co 1 f.. 
nt), 



(2.14) 



Equation (I2.14p . an important content of the present 
model, is independent of S!{oj). This expresses how the 
dissipative kernel 7(<) depends on the response function 
ip{t) of the medium due to the external noise e(t). 

If we assume that t(t) is a J-correlated noise, i.e., 
{e{t)e{t'))e = 2DS{t — t') then the correlation function 
of TT{t) will be 

{n{t)n{t')) = expHt - <'|/tc) (2.15) 

where we have neglected the transient terms (t, t' > Tc). 
This equation shows how the heat bath dresses the ex- 
ternal noise. Though the external noise is a 5-correlatcd 
noise, the system encounters it as an Ornstein-Uhlenbcck 
noise with same correlation time of the heat bath but 
with an intensity depending on the coupling kq and the 
external noise strength D. On the other hand, if the 
external noise is an Ornstein-Uhlenbeck process with 
{e{t)e{t'))e = {D/T')exp{-\t - t'l/r'), the correlation 
function of Tr{t) is found to be 



Dclnl r' 



{r'/Tcf 



-exp 

r' 



\t^t'\ 



1 

— exp 



\t-t'\ 



(2.16) 



where we have neglected the transient terms. If the ex- 
ternal noise-correlation time be much larger than the in- 
ternal noise-correlation time, i.e., t' ^ Tc, which is more 
realistic, then the dressed noise is dominated by the ex- 
ternal noise and we have from (j2.16p 



exp 



\t-t'\ 



(2.17) 



On the other hand, when the external noise correlation 
time is smaller than the internal one, we recover (|2.15[) . 



III. GENERALIZED FOKKER-PLANCK 
DESCRIPTION AND KRAMERS' ESCAPE RATE 

To start with we consider the internal dissipation is 
Markovian (i.e. Tc and the internal noise is Gaussian 
(5-correlated) 



7(t) = 2f3S{t - t') where /3 : 



(3.1) 



Consequently the generalized Langevin equation 
duces to 



dV{x) 

dx 
dV{x) 

dx 



(3{g'{x)fv + g'{x){f{t)+iT{t)} 
T{x)v + g'{x){f{t)+n{t)} (3.2) 



where 



nx)^p{9'{x)f 



(3.3) 



Using van Kampen's cumulant expansion method [29| . 
the Fokker-Planck equation, valid for small correlation 
time, corresponding to the above generalized Langevin 
equation is obtained as (see Appendix A) 



dP{x,v,t) 
di 



dP 
dx 

-[T{x)v + V'{x) - 2g'{x)g"{x)J, 



+A 



B 



d^p 

dvdx 



T{x)P 



dP 
(3.4) 



where 



A = {g'{x)fle - T{x){g'{x)fj, and B = {g'{x)fje, 

(3.5) 

and Ip and J„ are defined as 



{£,{mt-r))dT, 



Je = / T{£.{mt~T))dT. 
'0 



(3.6a) 
(3.6b) 



In equations (|3.6all3.6b|l . ^(t) is the effective noise term 
K(0 = + defined earlier. In deriving (|3.4p we 

have assumed that f{t) and e{t) are uncorrelated as they 
have different origin. Equation ^3.4^ is the first key result 
of this paper. It should be noted that when the noise 
is purely internal and the system-reservoir coupling is 
linear, equation l|3.4p reduces to the generalized Kramers 
equation Q (valid for small correlation time) . 

In Kramers' original treatment the dynamics of the 
Brownian particle was governed by Markovian random 
processes. Since the work of Kramers a number of au- 
thors [13, [m, m, [3^ have extended Kramers' analysis 
for the non-Markovian case to derive the expression for 
generalized escape rate. In order to allow ourselves a 
comparison with Fokker-Planck equation of other forms 
[31, 32, 34], we note that the diffusion coefficient in equa- 
tion (|3.4p is coordinate dependent. It is customary to get 
rid of this coordinate dependence by approximating the 
coefficients at the barrier top or potential well where we 
need the steady state solution of equation (|3.4p . One may 
also use mean field solution of equation (|3.4p obtained by 
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neglecting the fluctuation terms and putting appropriate 
stationary condition in the diffusion coefficient. The drift 
term in equation (|3.4p refers to the presence of a dressed 
potential of the form 



Rix) = V{x) ~ {g\x)fj. 



(3.7) 



The modification of the potential is essentially due to 
nonlinear coupling of the system to the nonequilibrium 
modes. Je is a non-Markovian small contribution and 
therefore the second term of the above equation may be 
neglected for small correlation time. For the rest of the 
treatment we use Tl(x) ~ ^ix)- For a harmonic oscillator 
with frequency wq, = uJqx'^ /2; the linearized version 
of Fokker-Planck equation is represented as 



dP 
'dt 



dP 



dP 



dv 



where 



dvdx 
(3.8) 



Ao = (5'(0))'/e - r(0)(g'(0))Ve and B„ = {g'{Q)fJe 

(3.9) 

are calculated at the bottom of the potential {x — 0). 
From equation (|3.9p we have 



(3.10) 



Ao = {g'{Q)fle - r(0)Bo. 



The general steady state solution of equation (|3.8|) be- 
comes 



Pst{x,v) 



■ exp 



,2^2 



2Do 2{Dq + Bo 



(3.11) 



where 



Ao 

r(o) 



(3.12) 



and Z is the normalization constant. The solution ()3.11|) 
can be verified by direct substitution in the steady state 
{dP{x = 0, v) /dt = 0) version of the Fokker-Planck equa- 
tion (|3.4p . namely 



-V— 1- TPst + [Fw + ujqx\— h Ao 



+5( 



dx 
dvdx 



dv 



dv"^ 



= 



(3.13) 



The distribution (j3.1ip is not an equilibrium distri- 
bution. In absence of the external noise e{t) it re- 
duces to the standard thermal Boltzmann distribution, 
exp[— (i;^ -|- V{x))/kBT]. Thus the steady state distribu- 
tion for the nonequilibrium open system plays the role 



of an equilibrium distribution of the closed system which 
may however be recovered in the absence of the external 
noise. 

We now turn to the problem of decay of a metastable 
state. In Kramers' approach [3], the particle coordinate x 
corresponds to the reaction coordinate, and its values at 
the minima of the potential well V{x) separated by a po- 
tential barrier, denotes the reactant and product states. 

Linearizing the motion around the barrier top at a; = 
Xb, the steady state {dP{x = Xb,v)/dt — 0) Fokker- 
Planck equation corresponding to equation p.4p reads 



dPst 
dy 

d^Pst 
dvdy 



2 dP,t 
= 0, 



r{xb)-^ivPst) + A, 

dv 



d^p, 



St 



d^P where 



y = x-Xb, V{y) = Eb 



^ 2 2 



> 



dv'^ 
(3.14) 

(3.15) 



and the suffix '6' indicates that all the coefficients are 
to be calculated using the general definition of A and B 
p.Sp at the barrier top. It is interesting to note that for 
linear coupling, we can extend our analysis for arbitrary 
correlation time and in such a case, the barrier dynamics 
would have been governed by the Fokker-Planck equation 
of Adelman's form (3^ . 

To derive a nonvanishing diffusion current across the 
barrier top Kramers f3\ considered Pb{x, v) to be the equi- 
librium Boltzmann distribution {exp[—{v'^+V{x))/kBT]) 
multiplied by a propagator F{x, v) and used it to solve 
the Fokker-Planck equation. In our model the equilib- 
rium distribution should be replaced by the steady state 
distribution which depends on the local nature of the po- 
tential at the barrier top and the effective temperature 
like quantity for the nonequilibrium open system. How- 
ever, in the absence of the external noise e(t) the steady 
state distribution reduces to the equilibrium Boltzmann 
distribution. Following Kramers 0] , we thus assume that 
the nonequilibrium steady state probability Pb{x, v) gen- 
erating a nonvanishing diffusion current across the barrier 
is given by 



Pb{x,v) 
where 



exp 



V 

2Db 



Dh = 



Vix) 



Db + Bb 



r{xb)' 



F{x,v), (3.16) 



(3.17) 



with 



Vix) 



Eq + -^ujqX , near the bottom 
Eb + -uJbi^ ~ , near the top. 
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The expression (|3.16l) denoting the steady state distribu- 
tion is motivated by the local analysis near the bottom 
and the top of the potential. For a stationary nonequilib- 
rium system, on the other hand, the relative population 
of the two regions, in general, depends on the global prop- 
erties of the potential leading to an additional factor in 
the rate expression. Because of the Kramers' type ansatz, 
which is valid for the local analysis, such a consideration 
is outside the scope of the present treatment. 

Following Kramers' original reasoning Q we then de- 
rive the barrier crossing rate K, for moderate to large 
friction regime (see Appendix B for detailed calculation) 



K 



A 



1/2 



exp 



(3.18) 

where Ej, is the barrier height of the potential and the 
parameter A is given by 



A = 



A 



Ab + aBb 



with 



r(.T,) + «(i + |^ 



2{Db + Bb) 



^y-T{xb)- ^T^Xb)+^io^ 



The strength of the external noise and the damping func- 
tion are buried in the parameters Z)o, So, -Db, Bb and 
A. Equation p.lSp is the second key result of this pa- 
per. Here we note that {Db -\- Bb)/kB in the exponential 
factor of (|3.18[) defines a new effective temperature char- 
acteristic of the steady state of the nonequilibrium open 
system and an effective transmission factor is contained 
in the prefactor controlling the barrier crossing dynam- 
ics. As expected both are the functions of the external 
noise strength and the coupling of the noise to the bath 
modes. In the absence of external stochastic modulation, 
e{t) = 0, equation (|3.18p reduces to standard Kramers' 
result Q, namely 



^0 
i'KLJb 




1/2 



exp 



f -Eb 



which can be verified with the explicit forms of the pa- 
rameters Z?o, -Bo, Db^ Bb and A with Je 0, /g = pksT, 
g{x) = X and D or kq equals to zero. 

As the understanding of the theoretical aspects of 
nonequilibrium statistical mechanics became better a 
vast body of literature have emerged modifying Kramers' 
original approach which are well documented in the re- 
view by Hanggi, Talkner and Borkovec [J]. In the fol- 
lowing we briefly discuss one of the issues from the post 
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(a) n=0.0 




FIG. 1: (color online) Variation of escape rate K with dis- 
sipation constant /3 with (a) jj, — and (b) /i = 0.5 using 
the parameter set bi — 1.0, Eh = 5.0, fcsT — 0.1, kq = 5.0, 
De = 0.1 and Te — 0.01, where solid line and circle represents 
analytical and numerical rates, respectively. 



Kramers development of the reaction rate theory which is 
closely connected to our work. Though Kramers' original 
approach was restricted by the Markovian assumption, 
however, in certain situations the memory effect becomes 
important and the generalized Langevin equation with a 
memory kernel must be accounted for. To the best of our 
knowledge following the seminal work in this direction 
by Grote and Hynes ^3^, Hanggi and Mojtabai [3l'| and 
Carmeli and Nitzan [32] have extended Kramers' work 
for an arbitrary memory friction and have found that 
the rate can often be larger than one would obtain from 
Kramers' orig inal approach. The analysis of Hanggi and 
Mojtabai [3l| was based on the non-Markovian Fokker- 
Planck equation of Adelman [34] for a parabolic poten- 
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FIG. 2: (color online) Variation of analytical rate constant K 
with kq for different values of the strength of nonlinearity ^ 
with bi = 1.0, Et = 5.0, fcsT = 0.1, = 1.0, De = 1.0 and 
Te = 0.01. 
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tial and they have used essentially the same approach of 
Kramers. The generalized Fokker- Planck a ppr oach have 
also been adopted by Carmeli and Nitzan [32| to derive 
the expression for the steady state escape rate in the 
non-Markovian regime. All the above mentioned theoret- 
ical approaches have been first supplemented with a full 
stochastic simulation by Straub,Borkovec and Berne [s^ 
where the authors explicitly studied the dynamics within 
the framework of non-Markovian generalized Langevin 
equation. Incorporation of memory effect in the above 
mentioned work gets reflected in the rate expression 



^0 




cxp 



knT 



where, 7 and uJb are long time limit of the memory kernel 
7(t) and the renormalized frequency at the top of the po- 
tential LObif), respectively. In this paper we have extended 
the above mentioned approaches [3ll, [H, [11] for state de- 
pendent diffusion in the non-Markovian regime to obtain 
a generalized steady state escape rate when the bath is 
modulated by an external stochastic force. Though our 
treatment is valid for small correlation time, it incorpo- 
rates most of the characteristics of the non-Markovian 
state dependent diffusion process. 
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IV. SPECIFIC EXAMPLE: HEAT BATH 
DRIVEN BY EXTERNAL COLORED NOISE 

As a specific example we consider that the heat bath is 
modulated externally by a colored noise e{t) with noise 
correlation 



, , . , f^^ De ( \t~t'\ 

(6(t)6(0)e = — exp -J ^ 

To V To 



(4.1) 



where and Te are the strength and correlation time 
of the external noise respectively. In addition to that we 
consider the internal noise /(t) to be white (i.e. Tc — > 0). 
The effective Gaussian-Ornstein-Uhlenbeck noise £^{t) = 
f{t) + Ti{t) will have an intensity Dji and a correlation 
time Tji given by [lB| 



Dr = 



TR 



mm)dt, 



D 



R Ja 



tmm)dt. 



Following the above definitions and using equation (|2.9 
we have (see equation (|2.17p ). 



TR 



cliksT 
Dr 



Dr 



-To- 



D,kI), (4.2a) 
(4.2b) 



FIG. 3: (color online) Variation of rate constant K with 
strength of nonlinearity (a) for different values of external 
noise strength (calculated analytically using ksT = 0.1) 
and (b) for different values of temperature fcsT (calculated 
numerically using De ~ 1.0). The values of the other param- 
eters used are bi — 1.0, Et, = 5.0, (3 = 3.0, kq = 5.0 and 

Te = 0.01. 



Consequently the functions A{x) and B{x) in equation 
becomes 



A{x) = ig'{x)fDR~r{x){g'{x)fTRDR, (4.3a) 
Bix) = i9'ix))^TRDR. (4.3b) 

From these equations we may evaluate the various param- 
eters to obtain the generalized escape rate from equation 



V. RESULTS AND DISCUSSIONS 

To study the dynamics we consider a model cubic po- 
tential of the form V{x) = bix^~b2X^ where bi and 62 are 
the two constant parameters with 61, 62 > 0, so that the 
activation energy becomes Ej, = Ab\/21b\. The nonlin- 
ear coupling function is taken to be g{x) = x+ (l/2)/ia;^, 
/i being a constant implying the strength of nonlinear- 
ity of the coupling function. We then numerically solve 
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FIG. 4: (color online) Turnover of rate constant K (numeri- 
cal) with dissipation constant P for (a) different values of ^ 
(for fcsT = 0.1) and (b)difrerent values of ksT (for /i = 0.5). 
The values of other parameters used are hi = 1.0, Eb = 5.0, 

= 1.0, Ko = 5.0, Te = 1.0. 



the Langevin equation (|3.2p using second order stochastic 
Heun's algorithm. To ensure the stabihty of our simula- 
tion we have used a small time step At = 0.001, with 
Ai/r/? ^ 1. The numerical rate has been defined as the 
inverse of the mean first passage time [ssf. The mean 
first passage time has been calculated by averaging over 
10,000 trajectories. The value of other parameters used 
are given in the respective figure. 

One of the result of Kramers' theory is that K varies 
inversely in the intermediate to strong damping regime. 
In Fig. I we have plotted the rate constant K against 
the damping constant /3 = Cg in the moderate to large 
damping region where our theory is valid and we compare 
the theoretical result (|3.18p with the numerical simula- 
tion data for two different values of the nonlinear cou- 
pling parameter /i. It is observed that the agreement 
between the theoretical prediction and numerical simula- 
tion is quite satisfactory. In Fig. 2 we plot the variation of 
rate constant K, obtained from theoretical result (j3.18p . 
with external coupling constant Kq for various nonlinear- 
ity parameters. We observe that for a given /z, the rate 
increases nearly linearly and for a particular value of kq, 
increase in /i causes decrease in rate, which is also ob- 
served from the Fig. 3(a) (where we have plotted K vs. 
fj, for different values of De from analytical result) and 



Fig. 3(b) (where the same has been observed numerically 
for different temperatures). 

In his dynamical theory of chemical reaction, Kramers 
identified two distinct regimes of stationary nonequilib- 
rium states in terms of dissipation constant /3. The es- 
sential result of Kramers' theory is that the rate varies 
linearly in the weak damping regime and inversely in the 
intermediate to strong damping regime. In between the 
two regimes the rate constant as a function of damping 
constant exhibits a bell-shaped curve known as Kramers' 
turnover [ij. In the traditional system reservoir model 
the dissipation and fluctuations are connected through 
the fluctuation-dissipation relation. A typical signature 
of this relation can been seen through this turnover phe- 
nomenon. Whereas for a thermodynamically open sys- 
tem where the heat bath is modulated by external noise, 
both the dissipation /3 and response function ip depend 
on the properties of the reservoir. Due to this connection 
between the dissipation and external noise source equa- 
tion (|2.9p , plays the typical role of thermodynamic consis- 
tency relation, an analog of the fluctuation-dissipation re- 
lation for a thermodynamically closed system, for which 
one can expect turnover feature for this open system. In 
Fig. 4(a) and Fig. 4(b) we have plotted the rate constant 
K obtained from Langevin simulation for a wide range 
of damping constant for different values of nonlinearity 
parameter fi and temperature, respectively. The figures 
apart form demonstrating the turnover of the rate con- 
stant with the variation of the damping constant, shows 
a shifting of maxima towards left i.e., weak damping 
regime with the increase of /x and also is consistent with 
Fig.3. 

While observing the variation of the rate constant K as 
a function of correlation time (re) of the external noise for 
different (Fig. 5(a)), different temperature (Fig. 5(b)) 
and for different nonlinearity parameter /i (Fig. 5(c)) we 
find an interesting result. In all cases K passes through 
a maxima, then decreases and ultimately becomes inde- 
pendent of Te for large values of Te . In short the rate of 
barrier crossing exhibits a resonance behavior with the 
correlation time of external noise which is responsible for 
fluctuation of barrier height. The above resonance phe- 
nomenon is known as resonant activation (RA) j3^. So 
far RA has been observed due to the barrier fluctuation 
as a result of direct driving of correlated noise to the sys- 
tem. In this model the barrier fluctuation occurs (|3.2p 
due to the driving of nonlinearly coupled heat bath with 
the system by correlated noise. So the RA occurs as a 
result of correlated noise driven bath, which is interesting 
and new feature, instead of direct system driving. The 
immediate experimentally observable situation could be 
if we consider a simple unimolecular isomerization reac- 
tion in a photochemically active solvent under the in- 
fluence of fluctuating light intensity (see section I), the 
reaction rate can be enhanced by tuning the correlation 
time of the fluctuating light field. It is also interesting 
to note that our nonlinear coupling model which yields a 
state dependent diffusion may have an important conse- 
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FIG. 5: (color online) Numerical variation of rate constant K 
with correlation of external noise Te for (a) different values of 
strength of external noise De (for ksT — 0.1 and fi = 0.5), 
(b) different values of ksT (for = 1.0 and fi — 0.5) and (c) 
different values of fj, (for ksT — 1.0 and De = 1.0). The values 
of other parameters used are 6i = 1.0, Et = 5.0, /? = 5.0 (for 
(b) /3 = 3.0) and Kq — 5.0. Note the logarithmic abscissa in 
(a) and (b). 



quence in the generation of current for a Brownian par- 
ticle moving in a periodic potential without any external 
bias. Because of its extraordinary success in explaining 
experimental observations on biomolecular motors active 
in muscle contractions, the state dependent diffusion has 
attracted wide attentions in recent years [32| • 



response function of the medium due to external noise 
can be established. We then followed the dynamics of 
the open system in the spatially-diffusion limited regime 
and derived the Fokker-Planck equation (corresponding 
to the Langevin equation with space dependent dissipa- 
tion and multiplicative noise) with space dependent diffu- 
sion coefficient containing the effective temperature like 
quantity. We then derived the generalized Kramers' es- 
cape rate for moderate to strong damping regime. From 
the point of view of the realistic situation we considered 
the special case where the internal noise is white and 
the external noise is colored and have calculated the es- 
cape rate for a model cubic potential. We also numer- 
ically simulate the Langevin equation and observe that 
the theoretical prediction agrees reasonably well with the 
numerical result. Dependence of the rate K on various 
other parameters have been studied and it is observed 
that the enhancement of rate is possible by tuning the 
correlation time of external noise. The creation of a 
typical nonequilibrium open situation by modulating a 
bath with the help of an external noise is not an uncom- 
mon phenomena in applications and industrial process- 
ing. The external agency generating noise does work on 
the bath by stirring, pumping, agitating etc to which the 
system dissipates internally. In the present treatment we 
are concerned with a nonequilibrium steady state that 
signifies a constant throughput of energy. We believe 
that these considerations are likely to be important in 
other related issues in nonequilibrium open systems such 
as thermal ratchet and related problems. 

So far, in this paper, we have considered a linear cou- 
pling in the interaction between the heat bath and the 
external driving noise, Hi„t- It will be interesting to see 
how the dynamics changes when Hint is nonlinear. Stud- 
ies have been made recently on anomalous diffusion in 
presence of correlated external noise [38*1 . Our present 
methodology may be extended to investigate the trans- 
port process when the bath is modulated by two corre- 
lated external noise. In our future communication we 
would like to address such issues. 



VI. CONCLUSION 

Based on a system reservoir microscopic model where 
the system is nonlinearly coupled to a heat bath which 
is modulated by an external, stationary Gaussian noise 
with an arbitrary decaying correlation function, we have 
generalized the Kramers' theory and have numerically 
analyzed the model to calculate the steady state rate 
of escape from a metastable well. We have shown that 
since the reservoir is driven by the external noise and the 
dissipative properties of the system depend on the reser- 
voir, a simple connection between the dissipation and the 
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APPENDIX A: DERIVATION OF 
FOKKER-PLANCK EQUATION 

Equation ()3.2p can be written in the following form 

Ml = Gi{ui,U2,t;f(t),Tr{t)) (Ala) 
U2 = G2{ui,U2,t;fit),Tr{t)) (Alb) 

where we have used the following abbreviations 

ui — X (A2a) 
U2 ^ V (A2b) 

and 



comes after it with respect to u. Making use of one of the 
main results for the theory of linear equation of the form 
(|A6|) with multiplicative noise [2§|, we derive an average 
equation for p [{p) = P{u,t), the probability density of 
u{t)] 




/•OO 

Ao + a^ / dT{Ai{t) exp{TAo)Ai{t — t)) 
Jo 



xexp(-rAo)}P. (A7) 

The above result is based on second order cumulant ex- 
pansion and is valid for the rapid fluctuations with small 
strength where the correlation time Te is short but finite 
0, I.e., 



Gi ^ V 



dVjx) 
dx 



(A3a) 

-V{x)v + g'{x){f{t)+Ti{t)}{Kih) 



The vector u with components ui and U2 thus represents 
a point in a two-dimensional "phase space" and the equa- 
tion (jAlallAlb]) determines the velocity at each point in 
this phase space. The conservation of the points now as- 
serts the following linear equation of motion for density 
p{u, t) in phase space [2§| 



dUn 



or more compactly 



dp 
dt 



(A4) 



Our next task is to find out a differential equation 
whose average solution is given by (p) [2§| where the 
stochastic averages has to be performed over two noise 
processes f{t) and e{t). V • G can be partitioned into 
two parts: a constant part V • Go and a fluctuating part 
V • Gi(i), containing these noises. Thus we write 



V • G{u, t, f{t)Mt)) - V • Go(m) -f aV • Gi(u, t, fit)Mt)) 

(A5) 

where a is a parameter (we put it as an external pa- 
rameter to keep track of the order of the perturbation 
expansion in aTg, Tg is the correlation time of the e(t), 
we put a = 1 at the end of the calculation) and also note 
that {{Gi{t))) = 0. Equation ||X4| therefore takes the 
following form 



p(-u, t) = [Aa + aAi)p{u, t) 



(A6) 



where = - V • Gq and Ai = - V • Gi. The symbol V 
is used for the operator that differentiate everything that 



(Ai(i)Ai(i')) =Ofor \t-t'\> 



Equation (|A7[) is exact in the limit Tg trends to zero. 
Using the expansions for Aq and Ai we obtain 



dP 
'dt 



/■OO 

-W-Go + a"^ dT(V • Gi {t) exp(-TV • Gq) 
Jo 



; V • Gi (t - r)) exp(TV • Go)} P. 



(A8) 



The operator exp(— tV • Gq) in the above equation pro- 
vides the solution of the equation 



u,t) 



dt 



-V • Go^(M,t) 



(A9) 



signifies the unperturbed part of p) which can be 
found explicitly in terms of characteristic curves. The 
equation 



ii = Go(u) 



(AlO) 



for fixed t determines a mapping from u{t = 0) to u(t), 
i.e., u u'^ with the inverse (u'^)^'^ = u. The solution 
of equation (|A9|) is given by 



d(u-*) 



^(u,t) =^(u-*,0) 



d{u) 



exp(-tV • Go)^(u,0), 
(All) 

\d{u *)/(/(«) I being a Jacobian determinant. The effect 
of exp(— <V • Go) on $^(m) is given by 



exp(-iV • Go)^{u, 0) = ^(m-*, 0) 



d(u-*) 



diu) 



(A12) 



The above simplification when we put in equation (|A8 
yields 
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where 



dP 
'dt 



-Go + dr 



d{u-^) 



d{u) 

x(Gi(w,i)V_r • Gi{u-\t-T)) 
d{u) 



d{u-^) 



P. 



(A13) 



where V_t- denotes differentiation with respect to w_t-. 
We put a = 1 for the rest of the treatment. We now 
identify 



A = {g'{x)fle - T{x){g'{x)fj, and B = (.g'(x))Ve, 

(A18) 

and /„ and J„ are defined as 



le = / {mi{t~T))dT, (A19a) 

JO 
1*00 

Je = / Tmat-r))dT. (A19b) 
Ja 



Ui — X 

U2 = V 

Goi = V. Gn = (A14) 

G02 = -Tix)v-V'ix) 

G12 = g'{x){f{t)+n{t)}. 

In this notation equation (|A13p now reduces to 



f - -|:(.m|;{r(.). + r(.)}p 

x{g'(x--)(/(t - r) + nit - r))}])P (A15) 

where we have used the fact that the Jacobian obeys the 
equation [2§| 



^ 1 



d{x\v^) 



d{x, v) 



so that the Jacobian equals to e~'"*. 

As a next approximation we consider 
the "unperturbed" part of equation (|AlallAlb|l and 
take the variation of v during Tg into account to first 
order in Tp. Thus we have 



APPENDIX B: DERIVATION OF ESCAPE RATE 



Inserting (|3.16p in ()3.14|) . we obtain the equation for 
F{x,v) in the steady state in the neighborhood of Xb 



dF 

il + Bb/Db)v— 



Db 



Db + Bb 



UjI{x ~ Xb) 



dF d^F 
+ T{xb)v] — + A^TT + -^fa a a 



At this point we set 



u ^ V + a{x — Xb), 



- 0. (Bl) 



(B2) 



where a is a constant to be determined. With the help 
of transformation (jB2p , equation (jBip reduces to the fol- 
lowing form 



d^F 

{A, + aBb}^ 



Db 2 



Db + Bj 



uj-^{x ~ Xb) 



+ {T{xb) + a[ 1 + — I ;>w 



dF 
du 



0. (B3) 



Now let 



X ^—x — Tv, V ~ V + Ttv + tV'{x). (AI6) 
Neglecting terms ^(t^), equation (|A16[) yields 



d _ ^ ^ ^ ^ ^ 
dv~'^ dv dx 



Taking this into consideration equation (jAlSp can be 
simplified after some algebra to the following form 



dP{x,v,t) _ dP 
dt dx 

+ [T{x)v + V'ix) - 2g'{x)g"{x)J,]^ 



+A 



Q2p Q2p 

B 



dvdx 



T{x)P (A17) 



Db 



Db + Bb 



ujbix - Xb) + ■{ r{xb) + a[l + jj- 



(B4) 

where A is another constant to be determined later. By 
virtue of the relation (|B4p . equation (|B3p becomes 



d'^F ^ dF 
dv? du 



where 



A 



A 



Ab + aBb 



(B5) 



(B6) 



The two constants A and a must satisfy the simultaneous 
relations 
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-Xa 
-A 



D 



b 2 



This implies that the constant a must satisfy the 
quadratic equation 



Db + Bb 2 , J., X 



D. 



Db ' ' Db + Bb 

which aUows the sohitions for a as 



b 2 n 



Db 



a± 



-T{xb) ± ^T^{xb)+ALu^ 



2{Db + Bb) 
The general solution of equation (jB5|) is 



(B7) 



F{u) = F2 exp f —j dz + Fi, 



(B8) 



where Fi and F2 are constant of integration. We look 
for a solution which vanishes for large x. This condition 
is satisfied if the integration in jBSj remain finite for 
|u| — > +00. This implies that A > so that only a_ 
becomes relevant. Then the requirement Pb{x,v) — > 
for X — > +00 yields 



To determine the remaining constant F2 we proceed 
as follows. We first note that as a; ^ —00 the pre- 
exponential factor in equation (jBlOp reduces to the fol- 
lowing form 



F 



2[-- 



1/2 



(B12) 



We then obtain the reduced distribution function in a; as 



Pb{x -00) = 27ri^2 



1/2 



exp 



Db + Bb 



(B13) 

where we have used the definition for the reduced distri- 
bution as 



Fix) 



P{x, v)dv. 



Similarly we derive the reduced distribution in the left 
well around a; « using equation p. lip where the lin- 
earized potential is V{x) = uo'^x'^ 



1 



P.,(.) = -v/2^exp^2(Z?o + i?o) 



(B14) 



Thus we have 

F{u) = F2 
and correspondingly 



exp 



2 



dz 



(B9) 



with the normalization constant \/Z given by 



Z 27r^Z?o(-Do + Bo) 



The comparison of the distribution (jBlSp and (|B14p 
near a; ~ 0, i.e.. 



Pb{x,v) 



F2 



X exp 



TT 

2A ^0 

t.2 



Az2 

exp I — I dz 



Vix) 



2Db Db + Bb 



(BIO) 



The current across the barrier associated with the 
steady state distribution is given by 



vPb(x — Xb, v)dv 



which may be evaluated using (jBlOp and the linearized 
version of V{x), namely V{x) ~ Eb — (1/2)ijj^(x — XbY 
as 



i = F2 



27r 



K + d: 



1/2 



Db exp 



-Eb 



Db + Bb 



(Bll) 



Pst{xo) = Pb{Xo) 



gives 



F. 



A 



1/2 



DbJ 2ny/2Tr{DQ + Bq) 



(B15) 



Hence from equation (jBlip . the normalized current or the 
barrier crossing rate K , for moderate to large friction, is 
given by 



K 



Db 



A 



1/2 



2n{Do + Bny/^ \1 + ADb 
where Eb is the potential barrier height 



exp 



-Eb 

Db + Bb, 
(B16) 
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